Consistent thermodynamic derivative estimates for tabular equations of state

Numerical simulations of compressible fluid flows require an equation of state (EOS) to relate the thermodynamic variables of density, internal energy, temperature, and pressure. A valid EOS must satisfy the thermodynamic conditions of consistency (derivation from a free energy) and stability (positive sound speed squared). When phase transitions are significant, the EOS is complicated and can only be specified in a table. For tabular EOS's such as SESAME from Los Alamos National Laboratory, the consistency and stability conditions take the form of a differential equation relating the derivatives of pressure and energy as functions of temperature and density, along with positivity constraints. Typical software interfaces to such tables based on polynomial or rational interpolants compute derivatives of pressure and energy and may enforce the stability conditions, but do not enforce the consistency condition and its derivatives. We describe a new type of table interface based on a constrained local least squares regression technique. It is applied to several SESAME EOS's showing how the consistency condition can be satisfied to round-off while computing first and second derivatives with demonstrated second-order convergence. An improvement of 14 orders of magnitude over conventional derivatives is demonstrated, although the new method is apparently two orders of magnitude slower, due to the fact that every evaluation requires solving an 11-dimensional nonlinear system.Comment: 29 pages, 9 figures, 16 references, submitted to Phys Rev

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead

Integrated Radiation Transport and Nuclear Fuel Performance for Assembly-Level Simulations

The Advanced Multi-Physics (AMP) Nuclear Fuel Performance code (AMPFuel) is focused on predicting the temperature and strain within a nuclear fuel assembly to evaluate the performance and safety of existing and advanced nuclear fuel bundles within existing and advanced nuclear reactors. AMPFuel was extended to include an integrated nuclear fuel assembly capability for (one-way) coupled radiation transport and nuclear fuel assembly thermo-mechanics. This capability is the initial step toward incorporating an improved predictive nuclear fuel assembly modeling capability to accurately account for source-terms and boundary conditions of traditional (single-pin) nuclear fuel performance simulation, such as the neutron flux distribution, coolant conditions, and assembly mechanical stresses. A novel scheme is introduced for transferring the power distribution from the Scale/Denovo (Denovo) radiation transport code (structured, Cartesian mesh with smeared materials within each cell) to AMPFuel (unstructured, hexagonal mesh with a single material within each cell), allowing the use of a relatively coarse spatial mesh (10 million elements) for the radiation transport and a fine spatial mesh (3.3 billion elements) for thermo-mechanics with very little loss of accuracy. In addition, a new nuclear fuel-specific preconditioner was developed to account for the high aspect ratio of each fuel pin (12 feet axially, but 1 4 inches in diameter) with many individual fuel regions (pellets). With this novel capability, AMPFuel was used to model an entire 17 17 pressurized water reactor fuel assembly with many of the features resolved in three dimensions (for thermo-mechanics and/or neutronics), including the fuel, gap, and cladding of each of the 264 fuel pins; the 25 guide tubes; the top and bottom structural regions; and the upper and lower (neutron) reflector regions. The final, full assembly calculation was executed on Jaguar using 40,000 cores in under 10 hours to model over 162 billion degrees of freedom for 10 loading steps. The single radiation transport calculation required about 50% of the time required to solve the thermo-mechanics with a single loading step, which demonstrates that it is feasible to incorporate, in a single code, a high-fidelity radiation transport capability with a high-fidelity nuclear fuel thermo-mechanics capability and anticipate acceptable computational requirements. The results of the full assembly simulation clearly show the axial, radial, and azimuthal variation of the neutron flux, power, temperature, and deformation of the assembly, highlighting behavior that is neglected in traditional axisymmetric fuel performance codes that do not account for assembly features, such as guide tubes and control rods